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Abstract 


Tha boundary-layer equations for two-dimensional 1nv*.ompress1b1e 
flow are Integrated numerically for the flow over a flat plate and a 
Howarth body. Injection is Introduced either Impulsively or period- 
ically along a narrow strip. Results Indicate that Injection perpen- 
dicular to the wall Is transmitted Instantly across the boundary layer 
and has rather little effect on the velocity profile parallel to the 
wall. The effect Is a little more noticeable for flows with adverse 
pressure gradients. Injection parallel to the wall results In fuller 
velocity profiles. Parallel and oscillatory Injection appears to In- 
fluence the mean. The amplitude of oscillation decreases with distance 
from the Injection strip but further downstream it Increases again In a 
manner reminiscent of an unstable process. 
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Introduction 

Blowing and suction has been an effective method of boundary-layer 
control for many years. One of the Important problems, namely how blow- 
ing and suction specifically affect the stability of a boundary layer 
has been addressed extensively. In this short report the effects of un- 
steady Injection are examined. However, this study concentrates on the 
response of a laminar boundary layer to unsteady vectored Injection. 

The Interpretation of the results and their significance In the study of 
flow stability are not discussed here. Some qualitative conclusions 
based on the fullness of the velocity profile can perhaps be drawn, but 
the Issue of stability Is not even addressed here. 

The purpose of this pilot project was to demonstrate the 
feasibility of the Idea of stability control by Injection. To this end 
the laminar boundary-layer equations were integrated numerically. 

During the short period of this Investigation It became obvious that 
this Is a rather limited approach to the problem. 

Unsteady Injection perpendicular and parallel to the wall was 
tested numerically. The disturbance was Introduced Impulsively In the 
flow over a flat plate. The effect of Impulsively started injection for 
the flow over a body with adverse presure gradient was also examined. 
Finally, oscillating Injection parallel to the wall was calculated. 

This report also contains suggestions on how to attack the problem In a 
more effective way through the boundary layer equations, as well as some 
comments on how the eal flow will respond and what to anticipate If the 
full Navler-Stokes equations are employed. Because of the limitations 
of time for this project as well as the limitation of this short report, 
a lot of Information about the actual numerical work has been left 
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out. However, the method employed has been described adequately In 
earlier publications which are referenced here. This report should 
therefore be considered a continuation of earlier published material. 
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!• Numerical Results 


Consider the unsteady boundary-layer equations in their dimen- 
sionless form 

3E * “ by. ^ ‘ ^ ^ 

(1.2) 

Note that v and y are already stretched by the square root 
Reynolds number. Numerical integration Is implemented in a 
plane in terms of Gdrtler Independent variables 

of the 
transformed 

« - /S 

(1.3) 


(1.4) 

where Ug(x.t) is the outer edge velocity. New dependent variables are 
also introduced. 


(1.5) 

'> "If f’ '' 

(1.6) 


A finite-difference method for solving the system of equations has 
been developed and adequately described in earlier publications 
(Telionis, Tsahaliis & Merle 1973). It employs an implicit scheme of 
integration and marches first in the direction of the axial distance 
Then it increments time and repeats the process by sweeping again the 
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physical space from the leading edge or a stagnation point. More 
details the reader w^ll find In earlier publications as referenced In 
Tellonis (1979). 

In this effort two body configurations Mere considered; a flat 
plate and a Howarth body. The second Is a fictitious body with an outer 
floM linearly decelerating according to the formula 

Ug » 1 - Ax (1.7) 

It Is equivalent to the flow over an ever steepening slope and therefore 
an Increasing adverse pressure gradient. 

On the flat plate blowing was Introduced In the domain 0.38 < IS < 
0.59, after a dimensionless time t > 0.5 had elapsed. The calculated 
vertical component at t ■ 0.48 Is plotted In Fig. 1.1 for different 
times. It appears that the entire boundary layer Is influenced Instant- 
ly and the vertical component continues growing with time. Moreover, 
the effect of blowing becomes stronger away from the wall. In other 
words, the disturbance actually grows with distance from the wall. For 
steady flow, this effect has been predicted by triple-deck theory 
(Nayfeh, Reed and Ragab, 1980). According to this theory the distur- 
bance of the v-component In the middle deck which is essentially what Is 
shown In Fig. 1.1 should be proportional to the mean profile 

» ■ - U(y) (1.8) 

where the function A 1$ a function of an appropriate axial scale and 
U(y) Is the Blaslus profile. Therefore, the disturbance on v grows with 
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distance y essentially like the u-component of the basic flow. 

The effect of blowing on the v-component further downstrem Is 
still evident as shown In Fig. 1.2, but not as violent as In Fig. 1.1. 
However, the u-component of velocity Is practically unaffected by blow- 
ing. This Is clear In Fij. 1.3. There Is a reasonable Justification 
for this fact which Is discussed In the following section. 

The effect of disturbances becomes more evident on boundary layers 
with advcr:~ pressure gradients, both favorable and unfavorable. In the 
present study the decelerating flow given by Eq. (1.7) with A « 0.12 was 
chosen which according to Howarth (1938) separates at approximately ^ * 
1, a fact that was verified by the present method. Separation of course 
Is not of Interest here but the sequence of profiles that lead to separ- 
ation are typical of a boundary layer developing against an adverse 
pressure gradient. A family of such profiles calculated with our steady 
flow program is shown In Fig. 1.4 for different ^-locations. The fuller 
profiles correspond to stations closer to the leading edge. 

Blowing perpendicular to the wall was Introduced again as before at 
a rate of \ ■ 0.4. The resulting u-components of the velocity are 
shown In Figs. 1.5 and 1.6 for times t ■ 0.25 and t ■ 0.75 
respectively. The effect Is again negligible. 

Unsteady blowing but parallel to the wall was also Introduced at 
the same stations and on the same flow field, namely the decelerating 
flow shown In Fig. 1.4. The results are shown In Fig. 1.7 and 1.8 for t 
■ 0.25 and t « 0.75 respectively. Again with u^^ ■ 0.05, the velocity 
profile is displaced from the origin, but little qualitative difference 
can be found. 

The method of plotting of Figs. 1.4 - 1.8 is perhaps not the most 
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appropriate because to compare profiles developing In time requires ob- 
servation of more than one figure. To bring out clearly the unsteady 
effect, results ««ere also plotted for a fixed station but at different 
dimensionless times. Blowing In the Interval 0.48 < ( < 0.S6 at a rate 
of > 0.2» Injected Impulsively after t • 0.25, resulted In the 
sequence of profiles shown In Fig. 1.9. The first three profiles 
correspond to i • 0.18, 0.28 and 0.38 and therefore undisturbed flow for 
all times, because blowing Is affected before ; > 0.48. A slight 
change with time In these profiles Is due to small adjustments In the 
solution, as the fictitious stagnation profile Is Iterated at each time 
step. The fam-llles of curves In 1.9d-h correspond to the undisturbed 
flow and the first three subsequent time steps with At « 0.25 
respectively. The profiles become less and less full as time 
progresses. In fact separ>ac1on is moving upstream and liito the domain 
of Integration and as a re>suit the solution does not converge at ( > 
0.88 and t ■ 0.75. This is vrtiy only three profiles are shown in Fig. 
1.9h. The case of suction has also been tested but the results are very 
similar with the opposite trend and are not shown to save space. 

Unsteady Impulsive Injection parallel to the wall at the same blow- 
ing strip and time and with a rate of u^/U^ « 0.05 results In the pro- 
files of Figs. l.lOa-b. The results for u^/Ug ■ 0.075 Indicate more 
clearly the development of the profiles as shown In Fig. 1.11. In these 
and all subsequent figures the profiles that correspond to stations ( ■ 
0.18, 0.28 and 0.38 are not shown because at these stations the flow Is 
not disturbed and therefore the profiles are Identical to those plotted 
In Fig. 1.9a, b and c. 

Finally the effect of oscillatory Injection parallel to the wall 


was tested with blowing through the same strip, namely 0.48 < ( < 0,56 
but with a sinusoidal magnitude, 

“w 

TT“ ■ B sin (lit 

e 

The dimensionless angular velocity (d was set equal to and 6 was given 
the values of 0.05 and 0.075. The results for these two values of B are 
shown in Figs. 1.12 and 1.13 respectively. The profiles plotted corre- 
spond to 16 times equally spaced In one period. They were obtained 
after three time periods have ellapsed to allow for a purely periodic 
motion to set In. Plotting In this maner results In overlapping of pro- 
files as can be seen In Fig. 1.13d. The results as plotted therefore 
can serve only as a qualitative guide or perhaps as Indicators of the 
envelope of the oscillating velocity profiles. 

It Is Interesting to note that further downstream from the distur- 
bance strip, the amplitude of oscillation quickly decreases, but as sep- 
aration Is approached. It once more Increases sharply. In fact even 
more Interesting Is the fact that near separation, the smaller blowing 
amplitude gives rise to a much more violent oscillation ':s shown In Fig. 
1.12f compared to Fig. 1.13d. 

What should be very Interesting to investigate here Is the mean of 
the oscillating profiles. Earlier studies (Tellonis and Romanluk, 1978) 
Indicate that periodic oscillations Imposed externally Induce a steady 
streaming component on the mean profile. Steady streaming (Riley, 1967) 
Is a net flow always In the direction of the adverse pressure gradi- 
ent. It could, therefore. Induce a fuller mean profile which could be 
more stable, provided that the frequency of the Injected oscillation Is 
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«May from the amplification range. 


2. Discussion 

The most characteristic property of Impulsive motions In fluids Is 
dwarfing of the convective terms compared to acceleration and diffu- 
sion. Blaslus (1908) recognized, for example, that for a boundary layer 
started Impulsively, the boundary- layer equations 


^ jv 
bx by 


0 


( 2 . 1 ) 


bu 

*bf 



+ 


bu 

■»y 


t 



b^ 

by2 


( 2 . 2 ) 


can be approximated by 


bx by 
bt bx 



(2.3) 

(2.4) 


For classical boundary conditions (no-si ip, no-penetration, etc.) the 
solution Is essentially the complementary error function 


u(x,y,t) - U,(x) erfc (2.5) 

/2t 

where U^(x) Is the outer flow velocity distribution after the Impulse. 
The Idea has been expanded upon first by the classical by now paper of 
Goldstein and Rosehnead (1938) and then by a large number of subsequent 
publications (see review articles of Stuart (1971) and Tellonis 
(1979). The profile of Eq. (2.5) contains the "local effect" In time 
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which has bean anticipated all along In tha present analysis* namely for 
very small times* the effect of the Impulsive motion Is confined to vary 
small distances away from the wall. The situation Is depicted schem- 
atically In Fig* 2.1. 

If a change Is Introduced Impulsively on an already developed flow* 
the dlsti'ibance 1s governed again by equations similar o Eqs. (2*1) and 
(2.2) (Sears, 1949; Tellonis, 1981). Moreover, If Instead of the outer 
flow* It Is the wall conditions 'hat change Impulsively, tne character 
of the solution remains tt«e same and the boundary conditions are read- 
justed, A classical problem that falls ^n this category Is tlw Rayleigh 
problem: an infinite flat plate given impulsively a finite velocity u^ 
In a fluid otherwise at rest. Then the boundary condltons are 


u ■ 0 

at t < 0 

(c^6) 

u ■ 0 

asy^- 

(2.7) 

u - u^ 

at y • 0* t > 0 

(2.8) 


and the solution becomes 

u ■ u^ erf (2.9) 

/2t 

where is the speed of the wall* as shown schematically In Fig. 2.2. 

Unsteady Injection parallel to the wall In an established boundary 
layer can be calculated by virtue of Eqs. (2.3) and (2.4) which govern 
the disturbance of the flow. 

Impulsive Injection perpendicular to the wall Is a completely dif- 
ferent problem. First of all, the counterpart of Rayleigh's problem 
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dots not exist. If «n Infinite flat plate Is given Impulsively a 
velocity perpendicular to Its plane rather than parallel to It, a unl- 
fom floM Is generated Instantly, assuming of course that the fluid Is 
Incompressible. A schematic of the flOM field Is shown In Fig. 2.3a. 

It Is due to the viscous efects that In Fig. 2.2 disturbances diffuse 
away from the wall acrordlng to the erf function. A more similar prob- 
lem to the one under Investigation here Is perhaps Injection through a 
finite opening on the plate. This then becomes essentially the problem 
of an Impulsively started jet, shown schematically In Fig. 2.3b. Again, 
at the first Instant, the solution should be Inviscid and therefore a 
uniform jet should extend Instantly all the way to Infinity. Only In 
finite time the velocity profile readjusts to a jet profile by virtue of 
viscosity and perhaps may start spreading sidewise. 

If flow is injected perpendicular to the wall In an established 
boundary layer, then It should Instantly penetrate across the entire 
thickness of the layer. This should be the behavior of the actual flow 
and this appears to be the response of the unsteady boundary-layer equa- 
tions, as discussed In the previous section. Since Inviscid motion In 
the y-directlon Is correctly modeled In the boundary layer equation, 
this behavior should have bec-n expected. However, for later times the 
model is not appropriate. The* diffusion of the v-component of velocity, 
namely the terms a^v/ax^ + are altogether missing from the 

boundary layer equations, Eqs. (2.1) and (2.2). This is because the 
second momentun equation has been eliminated altogether In the boundary- 
layer approximation. The reader should be reminded that It is due to 
diffusion terms that the erf or erfc functions appear in the impulsive 
motions given by Eqs. (2.5) annd (2.9). The effect of localized 
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disturbances In space and time Is therefore not present in the solution 
of the boundary-layer equations, but It should not be a feature of the 
real flow either. This Is because the establishment of a uniform 
unsteady Jet Is a purely Inviscid phenomenon, dominated for small times 
by continuity. A t • 0 It Is simply plug flow with speed throughout 
and a width equal to the width of the Injection opening. 

Steady blowing perpendicular to the wall cannot be analyzed via 
the boundary-layer equations unless the rate of blowing Is very small. 
If small blowing Is Introduced Impulsively, then there Is no justif- 
ication again for neglecting the second momentun equation since the 
term 8v/5t may be large, even If v, 9!/ax, etc., are small. A more 
correct set of equations then should be 


ax ay 


= 0 


( 2 . 10 ) 


du ^ 


„ au 


a^u 


M ^ 

ax 

ay ” 

"ax 


av 

■ s 


a^ 



at 

ay 

2 





ay 




( 2 . 11 ) 

( 2 . 12 ) 


The disturbance affects the pressure distribution In the vicinity of 
blowing. Pressure Is no longer constant across the boundary layer and 
should be considered a dependent variable. This set of equations can be 
Integrated numerically by the numerical method employed In the previous 
section. To this end, a second parabolic equation and a new dependent 
varlible will have to be added to the system. Inviscid boundary con- 
ditions can be Imposed at the edge of the boundary layer. The v- 
component of the velocity will not be matched with Its potential flow 
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counterpart as Is usually the practice In uninteracted calculations. A 
more appropriate procedure Is perhaps to consider the feedback from the 
disturbed potential floM. Allowing the boundary layer to Interact with 
the outer potential flow will make the nuiMrIcal solution equivalent to 
a Navler-Stokes solution or a tr1p1e>deck asymptotic analysis. 
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The v-component of velocity at C = 0.95 for blowing perpendicular to 
the flat plate wall. 
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The u-component of velocity at C = 0.48 for blowing perpendicular to 
the flat plate wall. 
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Fig. 1.4 Velocity profiles over the Howarth body given by Eq. 1.7 
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Fig. 1.5 The Influence of blowing 
pe’"pend1oular to the 
wall ; Time t » 0.25 


Fig. 1.6 The Influence of blowing 
perpendicular to the wall; 
Time t • 0.75 
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Fig. 1.9 Howarth flow with impulsiva injection perpendicular to the wall 
for time t ■ Q, Q.2S, 0.50, 0.75. 
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Fig. 1.12 continued 
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Fig. 1.13 continued 
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Schematic profiles of an impulsively started boundary 
layer; < 1^2 ^ ^3 










